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Q ' Abstract 

O , We have investigated the three-color Ashkin- Teller model (SAT), 

on the Wheatstone bridge hierarchical lattice, by means of a Migdal- 
KadanofF renormalization group approach. We have obtained the ex- 
J^ . act recursion relations for the renormalized couplings, which have been 

\^ \ used to investigate the phase diagram and to study the corresponding 

OO ■ critical points. The phase diagram, represented in terms of the dual 

transmissivity vector, presents four magnetic phases and nine critical 
J;^ I points. We have also numerically calculated the correlation length 

f— ^ ' (vt) and crossover (0) critical exponents, which show that seven of 

the critical points are in the Potts model universality class {q = 2, 4 

e 8). The remaining critical points are in a universality class which 

k>«( \ may belong to the Baxter's line. Our results are exact on the hierar- 

j_j ■ chical lattice used in the present work and the phase diagram can be 



considered as an approximation to more realistic Bravais lattices. 
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1 Introduction 

The N-color Ashkin- Teller (AT) model was introduced in order to investigate 
some aspects related to the non-universal critical behavior presented by the 
usual two-color AT model (2AT) [1]. Such non- universal behavior manifests 
itself in the two-dimensional 2AT model which exhibits a line of varying crit- 
ical exponents [2]. However, it was found that the non- universal behavior is 
a peculiar feature of the 2AT model. In Ref. pL] a series of techniques, such as 
first-order perturbation expansion, mean-field theory and Monte Carlo sim- 
ulations, were used to investigate the three-color AT (SAT) phase diagram 
in two and three dimensions. Since then, the SAT model has attracted some 
attention, mainly due to the richness of its phase diagram [21 IH El E]- In a 
recent study, the phase diagram of the SAT model was investigated by means 
of the exact Migdal-Kadanoff real space renormalization group treatment on 
a diamond hierarchical lattice [7] . In this work we perform a real space renor- 
malization group analysis of the SAT model on the Wheatstone hierarchical 
lattice [S] . Hierarchical lattices are constructed recursively in such a way that 
a large class of models can be solved exactly through a decimation procedure 
a la Migdal-Kadanoff real space renormalization group scheme [HIE]. 

In the present work we use a hierarchical lattice recursively constructed 
from an initial unit - the Wheatstone bridge - and a geometrical substitution 
rule as depicted in Fig. 1. Due to the lattice structure, we are restricted to 
consider only the ferromagnetic version of the SAT model under investiga- 
tion. Our results are exact in such a lattice and the phase diagrams may be 
considered as approximations to more realistic Bravais lattices. The plan of 
this article is as follows. In Section 2, we introduce the model and obtain 
the exact recursion relations for the renormalized couplings. In Section S we 
analyse the structure of the fixed points and determine the critical behavior 
of the model and the corresponding phase diagram. Finally, in Section 4 we 
summarize our main findings. 

2 The model and the renormalization group 
relations 

Let us consider the SAT model consisting of three Ising systems coupled by 
both four- and six-spin interaction terms. The Hamiltonian is given by 



/S'H = - ^ [i^2(o-j<Tj + 6i6j + TiTj) + K4{aicrj6i6j + atajnTj + OiOjTiTj) 

+KQ{(Ti(TjeiOjTiTj)], (1) 

where /3 = l/ksT, the sum is over the nearest-neighbor pairs of sites i and 
j and each spin variable {a, 9 and r) assumes one of two possible values 
independently: a = ±1, 9 = ±1 and r = ±1. Also, K2 = J2/kBT, K4 = 
J^/ksT and Kq = jQ/ksT represent dimensionless coupling constants. 

The SAT model described above will be studied on the hierarchical lattice 
generated as shown in Fig. 1. For this purpose we will employ a real space 
renormalization group approach which is known to be exact for hierarchical 
lattices and to give a good qualitative approximation for phase diagrams in 
real Bravais lattices [8] |9j. To treat the Hamiltonian of Eq. [1] we apply 
the renormalization group transformation indicated in Fig. 2. The renor- 
malization group recursive relations are obtained in two basic steps: (i) by 
the decimation of the spin variables associated with the intermediate sites 3 
and 4; (ii) followed by the requirement that the partition function should be 
preserved. 

As we have already mentioned in the Introduction, this kind of lattice 
becomes very suitable to investigate ferromagnetic order. For instance, in 
the case of the Potts model, it yields the exact transition temperature for the 
square lattice and good approximations to the critical exponents wherever it 
applies. Therefore, we will consider only stable ferromagnetic phases. This 
restriction impose some constraints that must be obeyed by the coupling 
constants as inequalities: 

K2 + Ki> 0, 3K2 + Kq>0 and K2 + 2Ki + Kq>0. (2) 

From now on we will make use of the transmissivity vector and its dual 
[lOj, as it was done in a recent work [7]. These vectors are very useful since 
they simplify some intermediate calculations and become very convenient 
to represent the phase diagram. The transmissivity vector is defined by its 
components [TD] 



I ^ ^{~2K2-4K4-2Ke) _ g(-4A-2 -4/4-4) _ ^{-(iK2-2Ke,) 
1 + 3g(-27^2-4X4-2X6) _|_ 3g(-4/f2-4K4) + g(-6i^2-2i^6) ' 

1 _ (,i-2K2-4.K4-2K6) _ ^(-AK2-4.Ki) _^ ^{-iiK2-2Kfi) 
1 + 3g(-2/^2-4X4-2X6) _|_ 3g(-4/f2-4K4) + g(-6i^2-2i^6) ' 



(3) 
(4) 



and 

I _ 3g{-2K2-4X4-2A'6) _^ 3g(-4i^2-4ii'4) _ g{-6K2-2/C6) 
^6 = ]_ _^ 3g{-2K2-4A-4-2A'6) + 3e(-4ii'2-4ii'4) + g(-6ii-2-2i^6) ' ^ '' 

whereas the dual transinissivity vector has components given by 

l^ ^ 1 + ^2-^4-^6 ^ ^(-27^-2-4^-4-2/^6) (■5') 



(-4K2-4i^4) (n\ 
1 -|- OL2 -|- Ot4 -|- Lq 

and 

g(-6X2-2i^6)_ (g) 

1 + 3t2 + 3t4 + te ^ ' 
After some calculations along the lines described in Refs. [7] and [8], we 
obtain the following recursion relations 

h = ^. (9) 

f; = §. (10) 



1 + 3t2 + 3t4 + h 


l-t2 


- 14 + te 


l + 3t2 


+ 3t4 + h 


l-3t2 


+ 3t4 - h 



and 



where 



f; ^ f . (11) 



+ 6^4^ + 3t/ + 6t/ + 12t/t/ + Qi^% 

+ 3^2^ + 6^2^ + 1, (12) 

+ 0^2 ^4 tg + 2t2^4 ^6 "I" ^^2 '^4 + ^'^2 ^4 

+ 6t2^t4^ + 4i^% + t^^ + i^), (13) 



A'2 — 2(^2 ^4*^6 + 2t2i4^6 + ^2 '^G + ^^^2 ^4 ^6 
+ 4^2 t4t6 + 2t2 t4tg + t^ + 5t4 + t^ 



and 

ATg = 2(^6^ + tV + 3^2%% + 6t2t4% 

+ 6i2%ie + Qi^i^ + Qi^i^ + ^i^i^) . (15) 

The recursion relations given by Eqs. [9]fTT] completely determine the 
critical exponents as well as the phase diagram for the SAT model. For 
instance, their fixed points and the renormalization fiows allow us to obtain 
the phase diagram. On the other hand, the critical exponents are obtained 
from the linearization around the fixed points. This will be discussed in the 
next section. 

3 Results and the phase diagram 

Using Eqs. [QlfTTj we were able to numerically determine the existence of four 
distinct phases separated by two-dimensional critical surfaces. The phases 
are characterized as follows (see also Fig. 3) 

• paramagnetic (P) — )■ corresponding to the trivial fixed point ^2 = 
t4 = tg = 1 and in which {a) = {6) = (r) = 0, {aO) = {ar) = {Or) = 
and {aOr) = 0; 

• ferromagnetic (F) — )■ corresponding to the trivial fixed point t2 = 

t4 = te = and in which {a) ^ 0, {6) ^ 0, (r) y^ 0, (aO) y^ 0, {ar) ^ 0, 
{9t) ^ and {aOr) =^ 0; 

• intermediate-one (Fi) — t- corresponding to the trivial fixed point 
t4 = 1, t2 = te = and in which (a) = {9) = (r) = 0, (aO) = (ar) = 
(er) = and {aOr) ^ 0; 

• intermediate-two (F2) — t- corresponding to the trivial fixed point 
tg = 1, t2 = t4 = and in which (a) = (9) = (r) = 0, (aO) ^ 0, 
(ar) ^ 0, {9t) ^ and {aOr) = 0. 



Each phase is characterized by its corresponding stable fixed point which 
are summarized in Table 1. We have also found nine additional non-trivial 



Table 1: Location of the trivial fixed points and their respective magnetic 
phases. 

Trivial fixed point (t2, ^4, to) Magnetic phase 

F (0, 0, 0) Ferromagnetic 

Fi (0,1,0) Intermediate-one 

F2 (0, 0, 1) Intermediate-two 

P (1, 1, 1) Paramagnetic 



fixed points, identified as the critical ones. The linearization around a non- 
trivial fixed point allows us to determine the associated critical exponents 
as 

UT = lnb/\n\i (16) 

and 

= lnAi/lnA2. (17) 

Here Ai > A2 are the two largest eigenvalues of the stability matrix and b is 
the scaling factor of the Wheatstone bridge. The critical points are presented 
in Table 2 with their respective critical exponents. We also note that except 
for the last two points the exact location of the critical points as well as the 
exact values for the critical exponents are given in closed form. 

The fixed points Ji, I2 and I3 have exactly one relevant eigenvalue and 
they belong to the Ising universality class (two-state Potts model). The 
unstable fixed point Ji is located on the surface which separates the F and 
F2 phases. In particular, it lies on the line given by t2 = t4 = which 
renormalizes onto itself - such kind of line will be called self-dual. We note 
also that this self-dual line connects two attractors - F and F2 - which 
corresponds to the ferromagnetic and intermediate-two phases, respectively. 
Similarly, the unstable fixed point I2 is located on the boundary between 
the P and Fi phases and lies on the self-dual line characterized by ^2 = ^6 
and t4 = 1. This line connects the P and Fi attractors which represent, 
respectively, the paramagnetic and intermediate-one phases. Finally, I3 is 
located on the boundary between the ferromagnetic and the paramagnetic 
phases, and it is an unstable fixed point on the self-dual line t4 = t2^, i^ = 12 
which connects the attractors F and P. 



Table 2: Unstable fixed points and their respective eigenvalues and critical 
exponents. 



Fixed 
point 



^8 

p(l) 
-M 

p(2) 

p(3) 
-'4 

h 

h 

h 

f/2 



(i2, ^4, i(i) 



/ 2v/2-l 2V2-1 2V2-l -> 
I 7 1 7 ! 7 ; 

(0, i 0) 



f- - 1) 

V3' 3' / 



(Ai,A2,A3) 



/- 32V2-t 
'> 17 



40-V2 40-V2 \ 
34 ' 34 > 



;/r 



(i, 0, 0) 
(S, 0, 0) 
(i,2V2-l,iS^) 
(0,0,^/2-1) (2\/2- 1,0,0) 

(\/2 - 1, 1, ^/2 - 1) (2^/2 - 1, 0, 0) 

(\/2 - 1, 3 - 2^/2, 5\/2 - 7) (2^/2 - 1, 29 - 20\/2, 71 - 50\/2) 
(0.27362, 0.27362, 0.18567) (2.33597, 1.27297, 0.86747) 
(0.28801, 0.22580, 0.28801) (2.33597, 1.27297, 0.86747) 



ln(2) 

ln(2) 
Nf) 
ln(2) 

Ml) 

ln(2) 

ln(2) 
ln(2\/2-l) 

ln(2) 
ln(2\/2-l) 

ln(2) 
ln(2v^-l) 

0.81698 
0.81698 



ln(22^|^) 



1„( 



40- V^^ 



lii{2v'2-l) 



3.51527 
3.51527 



From Table 2 we also note that the three fixed points labeled PJ , P4 

and P4 correspond to the ferromagnetic phase of the 4-state Potts model. 

The unstable fixed point PJ is located on the boundary between the F 

and Fi phases and lies on the self-dual line t2 = te = which connects 

(2) 
the attractors F and Fi. The unstable fixed point P4 is located on the 

boundary between the paramagnetic and the intermediate-two phases along 

the self-dual line given by t2 = ^4 and 1^ = 1 which connects the attractors P 

and F2. Finally, the unstable fixed point P4 is a multicritical point where all 

the four phases merge together and is located on the self-dual surface given 

by ^2 = U ■ ie, rightly at the intersection of two self-dual lines. One of such 

lines connects the critical points Ji and I2, while the other line connects the 

critical points PJ and PJ . Both P4 and P4 have one relevant eigenvalue, 

(3) 
whereas P4 has two relevant eigenvalues. 

The fixed point Pg belongs to the 8-state Potts universality class and lies 
on the intersection of two self-dual planes, t2 = ^4 and t2 = i^, i.e., it is 
located on the self-dual line 1 2 = t4 = tg which connects the attractors F and 
P. This fixed point is completely unstable and has three relevant eigenvalues, 
once Ai > A2 = A3, i.e., the lowest one is doubly degenerate. 

The nature of the remaining non-trivial fixed point f/i and U2 are quite 
distinct from the previously discussed. The stability analysis around these 
two fixed point reveals that they have the same set of three distinct eigen- 



values, two of which are relevant. Therefore, they have the same critical 
exponents vt and 0. These exponents do not correspond to any g-integer 
Potts universality class, as it has been verified in the case of the previously 
studied non-trivial fixed points. This peculiar behavior has been observed in 
a recent analysis of the SAT model on a diamond-like hierarchical lattice [?]. 
As far as their location is concerned, Ui lies on the merging line of the P, F 
and Fi phases, while f/2 lies on a line where the merging of the P, F and F2 
phases occurs. 

For the fixed points Ji, I2, I3, -P4 , -P4 , -P4 and Ps we have identified 
self-dual lines on which these points are located. This identification was 
crucial to determine the exact location for their respective coordinates as 
well as the associated critical exponents in closed form as presented in Table 
2. However, in the case of both Ui and f/2 we were not able to find out 
such identification of self-dual lines, and due to this fact we were limited to 
numerically determine the coordinates and critical exponents to these two 
points. In spite of this, a few words can be said about the points f/i and 
1/2- Note that the self-dual planes corresponding to ^2 = ^4 and ^2 = ^6 
shown in Figs. 4 and 5, respectively, are very similar to the findings reported 
in Ref. [13]. This resemblance may indicate that Ui and U2 belong to one 
of the infinity universality classes that exists along the Baxter's line. The 
motivation for this conjecture comes from the fact that for the 2AT model 
we can write the specific heat exponent a along the Baxter's line as 



where y is given by 



«-^^, (18) 

3-2|/' ^ ^ 



y = -cos-'{hexp{8K2)-l]}. (19) 

71 Z 

We must remark that K2 in Eq. (19) corresponds to the dimensionless four- 
spin coupling constant for the 2AT model, as defined in Ref. [15]. In the 
Baxter's line K2 is in the range [0, | In (3)], and as a consequence y is in the 
range [0, 1]. One can represent the parameter y along the Baxter's line as a 
function of the spatial dimension d and correlation length critical exponent 
ut (once duT = 2 — a), as follows y^ 

Sdur — 4 , , 



By substituting the value vt = 0.816 . . . , obtained for Ui and f/2 in the 
above expression we find for c? = 2 (corresponding to the square lattice) that 
y ~ 0.711... , whereas for df = In 5/ In 2 (the Wheatstone bridge fractal 
dimension) we get y ~ 0.942 .... Thus, in both d = 2 and df = In 5/ In 2 
cases the exponent ut for Ui and U2 may belong to the Baxter's line. How- 
ever, since we need two distinct critical exponents in order to determine the 
universality class we are left with a conjecture that should be probed by 
alternative methods. 

As a final remark, it is worth mentioning that all self-dual curves and 
surfaces we have found along the present work are also presented in a recent 
treatment of the SAT model on the diamond hierarchical lattice [7] . 

4 Conclusions 

In conclusion, we have studied the phase diagram and critical behavior of the 
SAT model. We have applied a Migdal-Kadanoff renormalization group tech- 
nique in a hierarchical lattice recursively constructed from the Wheatstone 
bridge. Due to the lattice structure we restrict ourselves to the ferromagnetic 
version of the model. Our results can be seen either as an approximation for 
the SAT on the square lattice, or as an exact result for the hierarchical lattice 
considered in the present work. We have found the presence of four magnetic 
phases, which are listed in Table 1 along with their attractors. We believe 
that the phase boundaries (the critical surfaces) obtained here represent a 
very good approximation to the exact results on the square lattice. We have 
also found nine non-trivial fixed points, which are listed in Table 2 along with 
their critical exponents. These nine non-trivial fixed points can be classified 
in four universality classes, three of them are associated to the Potts model 
with q = 2 states (Ji, I2 and I^), q = A states (PJ , P] and PJ ) and 
g = 8 states (Pg). Moreover, there are two non-trivial fixed points {Ui and 
f/2) which do not belong to the universality class of the Potts model. We 
should remark that all those results are qualitatively similar to what was 
found in a recent work (see Ref. [7j). More interesting, however, is the anal- 
ysis of the universality class associated to f/i and 1/2- Our numerical results 
show that the spatial dimension {d = 2 or df = In 5/ In 2) and correlation 
length critical exponent {ut = 0.81698) are compatible to one of the infinity 
universality classes that there are along the Baxter's line (see Eqs. 18-20). 
Therefore, we conjecture that both points Ui and U2 belong to one of the 



infinity universality classes that exist along the Baxter's line. 
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Figure 1: The graph used for the generation of the hierarchical lattice consid- 
ered in this work. The scaling factor 6 = 2 and the number of bonds m = 5 
furnish a fractal dimension df = In 5/ In 2. 
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Figure 2: Renormalization of the Wheatstone bridge hierarchical lattice. 
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Figure 3: Phase diagram of the SAT model in (^2,^4,^6) space. 
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Figure 4: Plane h = U of Fig. 3. 
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Figure 5: Plane ^2 = tg of Fig. 3. 
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